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We show that for superconducting Sr2Ru04 any unconventional pairing in the part of the Fermi 
surface with Ru 'idxy orbital character is weakly coupled to that with Ru {4da;2,4dyz} orbital 
character. This naturally gives rise to two disparate energy scales in the superconducting state which 
leads to novel low temperature properties in a variety of thermodynamic and transport properties 

1-^ ' and which would also account for the large residual density of states seen in specific heat and NQR 

rbr . measurements. 

t^^ ■ 74.70.-b,74.25.Bt,71.27.+d 

Sr2Ru04 provides the first example of a layered perovskite material that exhibits superconductivity without the 
presence of copper [|l|. Even though there is a close structural similarity with the high Tc materials, the electronic 
properties are very different. While it is clear that electron correlation effects are important in Sr2Ru04, the normal 



o 

^ • state near the superconducting transition is well described by a quasi-2D Landau Fermi liquid {e.g. the resistivity in 

.' [ all directions follows a T^ behavior for T ^ 50 K and the resistivity along and perpendicular to the c axis differ by a 

O I factor of 850 [|l|). Quantum oscillations show three Fermi surface sheets with a 2D topology that agrees well with band 

^ , structure calculations 1^. ft has been pointed out j|| that the mass enhancement is similar to that of '^He and that 

. ■ there is a metallic ferromagnetic phase in SrRuOs y^(the 3D analogue of Sr2Ru04). These observations indicate that 

j^ ' an odd-parity (/ = 1) superconducting state is likely g]. This is consistent with the lack of a Hebel-Slichter peak in 

Cj ^/Ti in NQR measurements [g] . A weak coupling analysis of the odd-parity state implies the gap should be of constant 

I magnitude |^. R is therefore surprising that specific heat |^ and NQR measurements Q reveal that approximately 

■ O 0.6 of the normal density of states remain in the superconducting phase in clean samples (those in which quantum 

oscillations were observed). As a consequence it has been proposed that an exotic non- unitary superconducting state 

similar to the "^He Ai phase is stabilized 03- In this scheme, the normal state quasiparticle energy spontaneously 

splits into two branches (one for spin up and one for spin down) upon entering the superconducting state. One of 

, • these branches is gapped while the other is not, leading to a residual density of states that is half the normal density 

Kj^ \ of states. 

QQ • Here we propose an alternative explanation for the large residual density of states. The electronic properties near 

^^ " the Fermi surface of Sr2Ru04 are determined by Wannicr functions with Ru dxy, d^z, and dyz orbital character P,fO|. 
^""^ ' We show that the quasi-2D nature of the electronic dispersion implies that the bands are derived from either the 
i.^ xy or the {xz, yz] Wannier functions and that the pair scattering amplitude between these two classes of bands 

, ' will be significantly smaller than the intraclass pair scattering amplitude for unconventional superconducting order 

^s^ \ parameters. R can therefore be expected that the gap on bands from one class is substantially smaller than that on 
"Y^ ■ bands from the other class. The presence of essentially gapless excitations for temperatures greater than the smaller 
Cd ' gap will appear as a residual density of states. Also, the two classes may favor different superconducting symmetries 
Q ' in which case a second superconducting transition will appear at low temperatures. 

Band structure calculations |^,|lO| give the correct shape of the Fermi surface, but predict an effective mass that is 
a factor of 4 smaller than that observed; indicating that strong coupling effects are important ||]. These calculations 
Q _ reveal that the density of states near the Fermi surface are due mainly to the four Ru Ad electrons in the t2g orbitals. 
O ■ There is a strong hybridization of these orbitals with the O 2p orbitals giving rise to antibonding tt* bands. The 
^ ' resulting bands have three quasi-2D Fermi surface sheets labeled a,l3, and 7 (see Ref. Q). The highly anisotropic 
nature of the Fermi liquid and the superconducting states suggests that the superconductivity essentially arises from 
intraplanar interactions, so we consider a single RUO4 plane. The Hamiltonian describing the band structure of a 
j^ . plane is 



ff= V i.,,.(R.-R,)4.,.,c,,,,-, (1) 



where Cn^i^s destroys an electron with spin s in the Wannicr function centered at R^ that transforms as the Ru v orbital 
{v = {xy, xz, yz}). Due to the Gz reflection symmetry about the center of the RUO4 plane ixy,a;z(R) = txy,yz(^) = 0. 
This implies that the 7 sheet of the Fermi surface can be attributed solely to the xy Wannier functions while the a 
and /3 sheets are due to a hybridization of the {xz, yz} Wannier functions. An effective Hamiltonian to describe the 
superconductivity is 



H =^ei{k)al^^^ai^i,^s+ ^ [V;,i'(k, k')aj'i^ ..al,;^ ,.,0;/ _k',s'ai',k',s + ^-c] (2) 

l,\i.,s /,r,k,k',s,s' 

where a^,k,s corresponds to the eigenoperators of Eq. (|^) and 

j,j',n,n' 

where t/(r, r') is an effective interaction and the spatial extent of the Wannier functions along the c axis restricts the 
integrations along z and z' to lie near the RUO4 plane. For the matrix elements Ky^^ and V-y^/j the symmetry of the 
Wannier functions under az can be exploited to write 

4C/(r^, z; r'^, z') = 2f/(r^, z; r^, z') - C/(r^, -z; r'^, z') - U{v^, z; r^, -z'). (4) 

The z dependence of the xy Wannier functions limits the integrations along the z direction in the V.y^a and V.y^p 
matrix elements to a distance on the order of 1/7 [ |ll| where / is the distance between two neighboring Ru ions. 
As a consequence, the lowest order term in a Taylor series expansion of Eq. ^ in z/|r^| and z'/|r'^l will give the 
largest contribution to V^ q, and V^^p for all but the onsite portion (Rj = R^/ = R,,„ = R„i') in Eq. 0. The lowest 
non-zero term is of second order in this expansion. Since the on-site contribution is independent of k and k' it 
does not contribute to the effective coupling constant for unconventional gap functions. It is therefore expected that 
the pair scattering amplitude between the 7 sheet and the {a, /?} sheets is significantly smaller that the intrasheet 
pair scattering amplitude (see Fig.P). Furthermore, since the Wannier functions forming the two classes of bands 
are of different symmetry, the intrasheet pair scattering amplitudes will in general be different. We assume that 
the superconducting state is odd-parity due to the considerations of Ref. Q. Note that the simplest tight binding 
approximation to the band structure (in which the Ru {dxzidyz} orbitals overlap only with neighboring O p — tt 
orbitals [pi) indicates that the gaps on the a and /3 sheets are the same magnitude for odd-parity pairing and we 
therefore assume that the gaps within this class have the same magnitude. 

We consider a model in which the three Fermi surface sheets have cylindrical symmetry and densities of states as 
in Ref. B. We use a weak coupling approach and in accordance with the above considerations take V;^;'(k, k') = 
Uijik ■ k' /{(kl}i{k'^)i>y^^ where {kf)i is the average of kf on sheet / and 

(^xy '^m ^m \ 
Um U U \ (5) 

Um U U ) 

where the matrix U operates on a basis with components that correspond to the Fermi surface sheets ^,a, (3 respec- 
tively. Introducing the gap matrix 

A,,^,, {I, k)^Y. ^.'' (1^' ^)Fs,.s, {I', k') (6) 

where i<"si,s2('j k) = (a(,k,siai,-k,s2): gives rise to a mean field Hamiltonian that is diagonal in the band index. For an 
odd-parity interaction the gap can be expressed as A(/, k) = i[d;(k) -crjcrj, jlj]. For unitary states (the case considered 
here) the quasiparticle excitations are given by Ei^i^ = (e^j^ + |d;(k)|^)^/^ and the gap equation is given by 

d.(k) = E^^^^%^^^^^tanh(/3i.^,,./2). (7) 

k ,/' 

Within weak coupling the transition temperature is Tc — 1.13ecexp[— l/Ama^;] where Xmax is the largest eigenvalue 
of the matrix with components Ui^i'{NiNiiy^^ and Ni is density of states of sheet I. It has been assumed that the 
cut-off frequency Cc is the same for all three bands. 

The superconducting order parameter is d;(k) — ^. ■ ci^ijkiXj/ {2kf)i which has a six fold degeneracy that is broken 
by spin-orbit coupling. The phases stabilized within weak coupling for the single band version of this model are the 
planar and the axial phases, both are degenerate within the approximations made above [n,f3|. Spin-orbit coupling 
will prefer one of these two phases and will fix the spin orientation of this phase to the crystallographic axes, leading 
to the classification in Ref. [B|. The quasiparticle excitation spectra for the possible phases are described by a gap 
of constant magnitude, so many properties will be correctly described by assuming that any one of these phases are 



stabilized. We assume that the Aiu phase, for which d;(k) — ci{xkx + yky)/[2k'^y/^, is stabilized. The resulting gap 
equation for the {q} then has the same form as that for isotropic superconductors generalized to include the presence 
of three bands |l|]. 

The interaction parameters Uxy, Um, and u remain to be specified. Earlier arguments imply Um ^ max{uxy, u) but 
the relation between Uxy and u remains unknown. Hund's Rule ferromagnetic correlations between the Ru dxz aud 
dyz orbitals may give rise to an increased odd-parity interaction for the {a, (3} Fermi surface sheets. Also, the 7 sheet 
is more 2D than the {a, (3} sheets, so fluctuations may lead to a greater reduction in the Tc for the 7 than for the 
{a,/3} sheets. These considerations indicate that u > Uxy, so for illustration purposes we consider this to be the case 
(though it cannot be ruled out that Uxy > u without a more detailed microscopic model). To show the qualitative 
behavior of this above model we take the density of states as measured in Ref. B] {Na : Np : N^ = 0.15 : 0.3 : 0.55) 
and the following values for the interaction matrix U; Uxy : Wm : u = 0.09 : 0.09 : 1.0 with uNp = 0.630. Using for the 
specific heat Ces 

and solving the gap equation yields the gaps and the specific heat shown in Fig.H. The presence of the small gap for the 
7 sheet gives rise to essentially gapless excitations for temperatures T > |d..y(k)| and this can give rise to the residual 
density of states observed experimentally. For temperatures below d^, this gap gives rise to the low temperature 
exponential decay of Ces/T to zero. Note that the density of states is split approximately evenly between the 7 sheet 
and the {a, /?} sheets. Consequently, the smaller gap lying in either the 7 or {a, /3} sheets gives good agreement with 
the magnitude of the residual density of states seen experimentally. To show how the smaller gap manifests itself in 
other properties we have calculated the London penetration depth and the thermal conductivity in the basal plane 
(shown in Fig. ). The London penetration depth is 
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which results from a simple extension of the standard BCS expression to include many bands. The thermal conduc- 
tivity in the single band case is derived in Ref. [ [l5[ and the suitable generalization to include many bands is 

«^-L(r) = -22^ -Tfrv^ij,—^ Ti^u (10) 

/,k ' 

with Tj^k = "nv.ilei.klZ-Ei.k where tj^ is the normal state relaxation time. This form is valid within the Born approxi- 
mation. It has been assumed that there is no interband scattering and that tn^i = tm- Note that rj^k does not have 
the same form as that for a conventional isotropic superconductor due to the odd-parity coherence factors |l6| , |l7| . 
In calculating these properties it has been assumed that the density of states corresponds to that of a clean system. 
However it may be the case that while the large gap will remain intact in the presence of impurities the smaller gap 
may be rendered gapless (though there will still be a coherent pairing amplitude on this Fermi surface sheet [[l2|). 

We have considered a model in which all the Fermi surface sheets favor the same superconducting symmetry. This 
model has two order parameters of the same symmetry (one (t/^i) for the 7 and one {'1P2) for the {a, /3} sheets) and can 
in principle have a second transition from a state in which (V'l, ^^2) = e'^(|ei|, ±|e2|) to a state in which time reversal 
symmetry is broken: (-01, ^2) = e*^(|ei|, e*'^|e2|) where 7^ 0, tt. An examination of the Ginzburg Landau coefficients 
found by a weak coupling analysis shows that the broken time reversal symmetry phase does not occur in this model. 
However, as was considered by Leggett for the two band conventional superconductor ||18| and more recently by Wu 
and Griffin in bilayer high Tc superconductors p9| , there will exist a collective excitation corresponding to fiuctuations 
into the broken time reversal symmetry phase (fluctuations of the relative phase of tpi and V'2). If all orbitals favor the 
same pairing symmetry then such a mode may appear below the single particle threshold. This mode is in addition to 
those that were predicted to exist due to the odd-parity symmetry in the presence of weak spin-orbit coupling [pj . It 
is also possible that due to the different symmetry properties of the the xy and the {xz, yz} Wannier functions the 7 
and {a, /3} sheets may favor different superconducting symmetries. In this case a second superconducting transition 
(as opposed to the crossover behavior shown in Fig. 0) is likely to occur due to the smallness of the pair scattering 
amplitude between these two classes of sheets. A low temperature broken time reversal symmetry phase is possible 
within this scheme |13]. 

In conclusion, we have presented a model for the superconducting transition in Sr2Ru04 in which the supercon- 
ductivity in the bands with Ru dxy orbital character and the bands with Ru {dxz, dyz} orbital character is weakly 



coupled. This model attributes the large observed residual density of states to thermal excitations across a secondary 
gap that is smaller than the primary gap driving the superconducting transition. This secondary gap should reveal 
itself in a wide variety of low temperature experiments on sufficiently clean samples. Also, within this model a sec- 
ond superconducting transition is possible. Experiments at very low temperatures are desirable to examine these 
possibilities. 
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FIG. 1. The vertex leading to the pair scattering ampUtude between the 7 sheet and the other two sheets of the 
Fermi surface. The effective interaction for any unconventional gap symmetry due to this vertex is small in relation 
to intrasheet interactions. 




FIG. 2. Specific heat, London penetration depth, and thermal conductivity as a function of temperature. The inset 
shows the magnitude of the gaps d-y and d^Q,.^} as a function of T /Tc- 



